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We report a first measurement of inclusive B — ^ X s rj decays, where X s is a charmless state with 
unit strangeness. The measurement is based on a pseudo-inclusive reconstruction technique and uses 
a sample of 657 x 10 6 BB pairs accumulated with the Belle detector at the KEKB e + e~ collider. For 



M Xs < 2.6 GeV/c 2 , we measure a branching fraction of (26.1±3.0(stat)l2;i(syst)l7;"(model)) x 10 

1-0.02 
-0.03 



and a direct CP asymmetry of Acp = — 0.13 ±0.04+°' {j 2 . Over half of the signal occurs in the range 



M Xa > 1.8 GeV/c 2 . 

PACS numbers: 13.25.Hw, 13.30.Eg, 14.40.Nd 



Decays of B mesons involving the b — ¥ s transition 
are an excellent tool for searches for physics beyond the 
Standard Model (SM). Theoretical treatments of these 
decays into exclusive hadronic final states, however, suf- 
fer from large uncertainties in the hadronization process. 
The uncertainties can be effectively reduced by leaving 
some of the final states in the calculation at the quark 
level, which corresponds to a measurement of an inclusive 
hadronic state X s of unit strangeness. 

Among such b — > s decays, those involving the 77 and 
7/ mesons exhibit unique properties due to interference 
between their underlying SU(3) octet and singlet com- 
ponents The CLEO collaboration reported the first 
measurement of inclusive B — > X s rj' with an unexpect- 
edly large branching fraction and an X s spectrum that 
peaks at high X s mass Q, a result confirmed in im- 
proved, higher-statistics measurements 0, H[ . Explana- 
tions included a large intrinsic cc component of the 77' 
the QCD anomaly mechanism Q that couples two 
gluons to the flavor singlet component of the 77', and 
also new physics sources Q. The first is disfavored by 
the lack of an enhancement of B — > rj c K relative to 
B — > J/tpK [8J, while the second is disfavored by a 
measurement of T(IS') — > r/'X Q, which indicates an 
rj'gg form factor that cannot explain the enhancement. 
A recent treatment [13] using soft collinear effective the- 
ory suggests that a measurement of the complementary 
process B —> X s r] can elucidate the possible contribu- 
tion from nonperturbative charm-penguin amplitudes or 
higher-order gluonic operators to both the 77 and 77' pro- 
cesses. CLEO performed the only previous search with 
an upper limit of B(B ->■ X s rj) < 4.4 x 10" 4 0. 



In this Letter, we report a measurement of B — >• X s r) 
using a sample of 657 x 10 6 BB pairs accumulated with 



the Belle detector at the KEKB e+e~ collider 11]. The 



Belle detector is a large solid-angle magnetic spectrome- 
ter and is described in detail elsewhere [l2j ]. 

We reconstruct candidate B mesons using a pseudo- 
inclusive method, with the X s composed of a K + or 
Kg(— > 7r + 7r _ ) and up to four pions, of which at most one 
is a 7r°(— > 77). This gives a total of 18 reconstructed 
channels and their charge-conjugates [3]. Candidate 
•q mesons are reconstructed in the 77 — > 77 mode from 
photons with E 1 > 200 MeV. The invariant mass of 
the 7-pair is required to lie between 520 MeV/c 2 and 
570 MeV/c 2 , or within 2a of the nominal mass. Wc 
veto an 77 candidate if either of its photons can be com- 
bined with another photon in the event to form a candi- 
date 7T°. To suppress background from radiative B de- 
cays, we require the energy asymmetry of the two pho- 



tons, defined as |i? 7l — E l2 \ 



\E lx + E J2 \, to be less than 



0.6. The momenta of 7/ candidates are recalculated us- 
ing the nominal 77 mass 14]. To suppress secondary 77- 
mesons from b — > c — > 77 chains, we retain only 77 can- 
didates whose center-of-mass (CM) momentum satisfies 
|p*| > 2.0 GeV/c. 

Charged pions and kaons are selected based on infor- 
mation from the time-of- flight, aerogel Cherenkov, and 
drift chamber dE/dx systems. Typical efficiencies to cor- 
rectly identify kaons (pions) are above 88% (98%), with 
misidentification rates for pions as kaons (kaons as pions) 
below 12% (4%). Kg candidates are required to have an 
invariant mass within 16 MeV/c 2 (4er) of the Kg mass 
and a displaced vertex from the interaction point. For 
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7T° candidates, each daughter photon is required to have 
energy greater than 50 (100) MeV in the barrel (end- 
cap) region and a shower shape consistent with a photon. 
The invariant mass of the photon pair must be within 
15 MeV/c 2 (2.5<r) of the ir° mass. The 7r° momentum 
is recalculated using the nominal 7r° mass. To suppress 
combinatorial backgrounds, we require 7r° candidates to 
have laboratory momenta greater than 300 MeV/c. 

Pions and kaons are combined to form an X s , and 
B meson candidates are formed from combinations of 
an r\ and an X s . A beam-constrained mass, M\> c = 
p* g \ 2 /c 2 and energy difference, AE = 
are calculated, where -Ebcam, p* B , and Eb 
are the beam energy, B momentum, and B energy, all in 
the CM frame. The signal is obtained using fits to Mb c 
with \AE\ < 0.1 GeV. 

We use a simulated signal Monte Carlo (MC) sam- 
ple [H consisting of B -> K-q for M x , < 0.6 GeV/c 2 , 
B -> K*r] for M Xs E [0.8,1.0] GeV/c 2 , and B -> X s r) 
in all other mass regions (Mx s € [0.6,0.8] GeV/c 2 , and 
M x , > 1.0 GeV/c 2 ). For the B -> A s ?7 component, 
fragmentation of the X s system into hadrons is simu- 
lated by PYTHIA [16(, assuming a model in which the 
X s mass spectrum is flat from the Kit threshold up to 
3.2 GeV/c 2 . We find an average of approximately nine 
B candidates per event, of which we select the candi- 
date with the lowest % 2 . This x 2 is defined as the sum 
of %a_b = {AE/<tae) 2 , where the resolution cta e is es- 
timated separately for each reconstructed mode, and, if 
available, a reduced-x 2 of a vertex fit that includes all 
X 8 daughter charged tracks except those used as part of 
a Kg candidate. After applying this procedure, we select 
the correctly reconstructed B in 56% of simulated events. 
The dominant background to B — > X a rj comes from 



continuum production of quark pairs, e + e 
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u,d,s,c). These events have a jet-like topology, and are 
suppressed relative to the spherical BB events using a 
Fisher discriminant 17[ formed from event shape vari- 
ables [H, [l9| . Further suppression is obtained by com- 
bining this Fisher discriminant with the cosine of the B 
flight direction in the CM frame and, when available, the 
displacement between the signal B and the other B in 
the event. This suppression is optimized as a function of 
b- flavor tag quality J2(j], and is approximately 34% effi- 
cient for the signal modes while suppressing over 99% of 
the continuum background. 

Decays of the type B — > X c rj and B — >• X c — > X s r], 
where X c is any state containing charm mesons, may 
have final states identical to the signal mode. We search 
among the candidate B decay products for combinations 
consistent with selected charm meson decays and veto 
the candidate if the mass of the reconstructed combina- 
tion is within ±2.5<r of the known mass. The modes 
and their veto widths are: D° — > Kmr^^n ), 13.5 
(44.5) MeV/c 2 ; D+ -4 Amr^Tr ), 12.5 (31.3) MeV/c 2 ; 
D" -> K s r), 31.3 MeV/c 2 ; L>+ ??7r+, 29.3 MeV/c 2 ; 



and r] c (lS) — > r)ir + Tr~ , 85.0 MeV/c 2 . We also veto events 
with an rj' —± ry7r + 7r _ candidate with an invariant mass 
M„^ within 100 MeV/c 2 of the nominal rj' mass. 

Signal yields are obtained using an extended unbinned 
maximum likelihood fit to M^ c in 200 MeV/c 2 bins of X s 
mass up to 2.6 GeV/c 2 . The probability density function 
(PDF) for the signal is taken as a Gaussian, with the 
mean and width determined from the appropriate signal 
MC sample (K-q, K*i], or X s rj) for the mass bin, with 
calibration factors taken from a B — > Dir control sample. 
All reconstructed modes are combined for the fit, and 
no attempt is made to separate correctly reconstructed 
B candidates and those with some missing or incorrectly 
attributed B daughters (self-cross-feed). Shapes for the 
charm contributions remaining after the vetoes are as- 
signed based on a MC sample of generic b — > c pro- 
cesses. Four separate PDFs are assigned for the largest 
charm backgrounds as identified in MC: B° — > D°r), 
B" -> D*°r), B" -> dW-tt+17, and B+ -> D^tt+t]. All 
other b — > c backgrounds are combined into another PDF. 
Each charm PDF consists of a Gaussian component to de- 
scribe the peaking in Mbo and an empirically determined 
parameterization (ARGUS function) [21| to describe non- 
peaking combinatorial contributions. The shape parame- 
ters are taken from the appropriate background MC sam- 
ple. Normalizations of the modes B° — > D^°rj are based 
on the previous Belle measurement [22]. The branch- 
ing fractions for the decays B — > D^nr] are unknown, 
so their normalization is determined by a simultaneous 
X 2 minimization based on the difference between the ex- 
pected and observed Mb c distribution of the events in 
all eight veto windows. The relative normalizations of 
D^tt + i] and D*~ir + T] are assumed to be the same, and 
likewise for the DW modes. The \ 2 technique is veri- 
fied by repeating the optimization over the B° — > D^°r] 
modes, for which the results are consistent with the previ- 
ous Belle measurement. This \ 2 1S also used to study sys- 
tematic errors on the normalizations of all charm PDFs. 
Normalization for the PDF that includes all other b — > c 
modes is fixed to the MC expectation. The remaining 
combinatorial qq backgrounds are modeled with an AR- 
GUS function. For the final fit, the signal yield and both 
the yield and shape parameter of the qq ARGUS PDF 
are allowed to vary. 

Rare B decay backgrounds are studied with a dedi- 
cated MC sample, and include contributions from B — > 
X a r)', B — > A s 7, and B — > Xdq. These expected yields 
are subtracted from the fit yield to give a final yield. 
The expected yields for B — > X s r\' and B — > X s ^ are 
based on the known branching fractions, and are found 
to be less than 0.5 events in each X s mass bin. The 
B + Tr + r) branching fraction is also known, and the ex- 
pectation is 5.2 events in the lowest bin of X s mass. We 
estimate the contribution from other B — > XdX\ modes 
by repeating the reconstruction and the fitting proce- 
dure but replacing the K + candidate of X s with a ir + 



candidate. Performing these fits on data and using a 
dedicated X^t] MC sample to estimate the rate to misre- 
construct Xd as X s , we estimate a total contamination 
of 19.1 ± 2.3 events from X^rj, distributed uniformly in 
the range M Xs G [0.6,2.6] GeV/c 2 . 

The fit to the full mass range, M Xa G [0.4, 2.6] GeV/c 2 , 
is shown in Fig. QJa) , and gives a background-subtracted 
yield of 1054 ± 54+]J. We also define a high mass re- 
gion, Mx 3 G [1.8,2.6] GeV/c 2 , where the summed yield 
is 233±34l]J. Significances are determined in each mass 
bin by convolving the likelihood function with a Gaus- 
sian of width determined by the systematic errors on the 
yield. The maximum likelihood, Anaxj an d the likelihood 
at a signal yield of zero, C$, are used to determine the 
significance, which is defined as — 21n(£ /£max)- The 
significance is 23 (7) for the full (high) X s mass range. 
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FIG. 1: (color online), (a) The Mb c distribution for the full 
mass range, Mx s £ [0.4, 2.6] GeV/c 2 . The points with errors 
correspond to the data, while the curves correspond to the 
overall fit PDF (solid red), the signal PDF (dashed magenta), 
the sum of all b — > c background PDFs (dotted green), and 
the combinatorial background PDF (dash-dotted blue), (b) 
Differential branching fraction, dB/dMx s , for B — > X s r). The 
three error bars correspond to statistical error only, statistical 
plus systematic error, and total error. 

Reconstruction efficiencies in bins of X s mass range 
from 6.5% to 0.1%, not including the branching frac- 
tion for rj — > 77; these results are based on the signal 
MC and assume equal production of B + B~ and B°B° 
at T(4S). Figure (Ub) shows the differential branching 
fraction as a function of Mj a . Table Q] gives the final 
results for each X s mass bin. For the full Mj a range, 
we sum the individual contributions and find the fol- 
lowing partial branching fraction B(B — > X s rj;Mx s 6 
[0.4,2.6] GeV/c 2 ) = (26.1 ± 3.0±^t£5) x 10~ 5 , where 
errors are statistical, (model-independent) systematic, 
and decay modeling. A large fraction of the inclu- 
sive signal occurs in the high mass region, where we 
find B(B -)■ X s <n;M Xa e [1.8,2.6] GeV/c 2 ) = (16.9 ± 
2.9 (state's (syst)^ (model)) x 10" 5 . 

The direct CP asymmetry is defined as Acp = (B~ — 
B+)/(B- +B+), where B+(B-) is the partial branch- 
ing fraction for B + or B° (B~ or B°). We measure 
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TABLE I: Measured background-subtracted signal yields 
(Ns), branching fractions (B), and CP asymmetry (Acp), 
for each Mx B range. Uncertainties on Ns are statistical. Un- 
certainties on B are statistical, systematic, and modeling, re- 
spectively. The uncertainties for Acp are statistical and sys- 
tematic. 



M Xs (GeV/c 2 ) 


N s 


23(10~ 6 ) 


-4cp(10~ 2 ) 


0.4- 


-0.6 


60 ± 12 


1.9 ± 0.4 ± 0.1 ± 0.0 


o f i -in i c\ 

—35 ± 18 ± 2 


0.6- 


-0.8 


15 ± 9 


0.9 ± 0.5 ± O.ll^o 


2 ± 40 ± 13 


0.8- 


-1.0 


250 ± 19 


17.0 ± 1.3 + ?o ± 0.0 


-4 ± 7 ± 2 


1.0- 


-1.2 


84 ± 14 


72± 1 2 + °i +0 A 


—26 ± 15 + 4 


1.2- 


-1.4 


146 ± 17 


15.8 ± 1.9 ± 1.0±i;S 


-22 ± 11±3 


1.4- 


-1.6 


137 ± 18 


20.8 ± 2.7±1;1±2.1 


-15 ± 12±3 


1.6- 


-1.8 


128 ±18 


28.2 ± 4.1 ± 2.1±|;f 


-25 ± 13±3 


1.8- 


-2.0 


64 ± 18 


24.4 ± 6.8±|;|±^ | 


-31 ± 26 ± 6 


2.0- 


-2.2 


86 ± 18 


42.4±9.lt^l^ 


34 ± 20± 3 


2.2- 


-2.4 


49 ± 18 


36.8 ± 13.5+^ +™ 5 


2 ± 32 ± 5 


2.4- 


-2.6 


35 ± 13 


65.1±23.4±? 2 5 8 ± 2 « 


-40 ± 36l:i 2 


0.4- 


-2.6 


1053 ± 54 


261 ± 30+^+^ 


-13±4l 2 


1.0- 


-2.6 


728 ± 48 


241 ± 30t^l 4 7 ? 


-15 ±6 ±3 


1.8- 


-2.6 


233 ± 34 


169 ± 29ii|ll 


± 14 ± 5 



this asymmetry in the subset of reconstructed modes 
in which the B flavor can be inferred from the final 
state. We correct the fitted CP asymmetry to account 
for events that are reconstructed with the wrong B fla- 
vor by multiplying the raw fitted asymmetry by a di- 
lution factor. This factor is estimated from the signal 
MC, and ranges from unity to 1.05. The bin-by-bin re- 
sults, as well as the results of separate fits for Acp over 
the full X s mass range and the range above the nar- 
row kaonic resonances (Mx s G [1.0,2.6] GeV/c 2 ), are 
shown in Table H For M x , G [0.4,2.6] GeV/c 2 , we find 
A C p = -0.13 ± 0.04t°;° 2 , with a significance of 2.6cr rel- 
ative to a null asymmetry. All Acp results that include 
the range Mx 3 G [0.4,0.6] GeV/c 2 are calculated with 
the assumption that the B + — > ir + i] backgrounds in this 
region contribute a CP asymmetry consistent with the 
existing measured world average [141 ] . 

Systematic errors on the fitted signal yield are domi- 
nated by PDF uncertainties. Uncertainties in the signal 
PDF parameters are studied using a B — > Dtt control 
sample, while those due to normalizations and shapes for 
the b — >■ c backgrounds are estimated using comparisons 
between the veto window \ 2 procedure and either MC 
expectations or, when available, previous measurements. 
Errors from the background subtractions are dominated 
by uncertainties in the estimate of backgrounds from 
B — > XdTj. Our estimates of these backgrounds may 
have included other small contributions, such as those 
from B — > X s ri, so we allow these estimates to vary by 
— 100%. Positive uncertainties are estimated from the 
difference in expected yields assuming a flat distribution 
of Xd events in X s mass versus those obtained from a 
MC study of cross-feed from Xd mass to X s mass. In 
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all cases the systematic uncertainties on the background- 
subtracted signal yields are at least a factor of two smaller 
than the statistical errors. 

The model-independent systematic error includes con- 
tributions from the signal yield, the selection efficiency, 
the number of BB pairs, and the r\ — > 77 branching 



A-cp by the measured uncertainty [14( (±100%) 



fraction [14|. For X s mass bins above 1.8 GeV/c , the 



errors on the signal yields from uncertainties in the PDF 
shapes (primarily for the charm PDFs) dominate with a 
contributed relative uncertainty of 7-18%. For the lower 
X s mass bins, the efficiency error is the largest contri- 
bution with a relative uncertainty of 5-6%. This error is 
the combination of individually determined contributions 
from control sample studies of the following: tracking, 
reconstruction of n°, 77, and Kg, particle identification, 
continuum suppression, and candidate selection. 

We define an additional error due to modeling of the X s 
system, which is studied in three parts. The first is due to 
the fraction of unreconstructed modes (e.g., modes with 
too many total or neutral pions or additional kaons) . We 
vary these fractions by ±30% of the PYTHIA expectation 
and use the differences in efficiency to estimate an Mx. 
bin-dependent uncertainty that rises with X s mass from 
zero to ±21.1%. The second is due to differences in the 
observed frequency of decay modes and those expected 
from PYTHIA. We find good agreement between data 
and MC in the relative amounts of charged and neutral 
B modes, modes with K s and those with K + , and modes 
with one or two total 7r's and those with three or four to- 
tal 7r's. However, we find a significant excess of modes 
without a 7T over those with a ir° , which we attribute to 
inaccuracies in the PYTHIA fragmentation. To quantify 
this uncertainty, we re-estimate the PYTHIA efficiencies 
with the fraction of ir° modes adjusted to match data, 
and use the difference between this value and the nom- 
inal efficiency to assign an error. This error is usually 
only negative, due to the higher reconstruction efficiency 
for modes without a 7r°, and is as large as —37% in the 
highest X s mass bin. The final component of the mod- 
eling uncertainty is due to the assumed X s mass spec- 
trum. We study the efficiencies of other Mx B signal MC 
samples where the spectrum rises toward high mass and 
assign errors based on the differences from the flat Mx s 
MC. Using these samples, we also study the fractions of 
self-cross-feed candidates that are reconstructed with an 
incorrect X s mass. These effects are small compared to 
the first two components of the modeling error. 

The systematic error on Acp includes contributions 
due to: uncertainties in the PDF parameters; possible 
detector and measurement biases, which are estimated 
from the measured Acp of the B — > Dir control sam- 
ple and the signal MC, respectively; uncertainty due to 
the signal model is studied by checking the fractions 
of events with incorrectly identified flavor using alter- 
native Mx s spectra; and possible contamination due to 
B — > ttt] (B — > XdTj) decays is estimated by varying their 



In summary, we report the first measurement of the 
inclusive process B — » X s rj, and find a partial branch- 
ing fraction of B(B X s r); M Xs € [0.4, 2.6] GeV/c 2 ) = 
(26.1 ±3.0(stat)±i'5 , (syst)±7;? (model)) x 10~ 5 . The mea- 
sured Mx 3 dependent branching fractions are consistent 
with the known B — > Kr\ and B — > K*(892)r) processes 
. In the high mass region, M Xa G [1.8,2.6] GeV/c 2 , 
which is above any significant contributions from pre- 
viously measured exclusive processes [24j |. we observe 
a signal with a 7a significance. We also measure the 
CP asymmetry of B —> X s r), both as a function of 
Mx s and for the full mass range, where we find Acp — 
-0.13 ±0.041^3- 

No theoretical prediction is currently 
available for the shape of the Mx s spectrum. However, 
the similarity in spectral shape to B — > X s rj' and the lack 
of strong suppression of the B —> X s rj branching fraction 
relative to the rj' mode imply that the origin of the large 
contribution in the rf mode is also common to the 77 mode 
[Toj ] , and disfavors rj specific mechanisms [H], Q . 
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